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ABSTRACT 

We proposed and solved a combined control malaria 
system of fifteen ordinary differential equations 
modeling the transmission dynamics of malaria 
between humans and mosquitoes. Since our aim is to 
minimize the number of exposed and infectious 
human, we discussed the disease free equilibrium and 
estimated the basic reproduction number using the 
next generation matrix method. The disease free 
equilibrium was asymptotically stable when the 
reproduction number is less than one and unstable 
when it is greater than one. Numerical results are 
provided using matlab software to confirm the 
analyzed results. Our findings were that malaria may 
be controlled using the combined control method, the 
insecticide treated bed nets, fumigation of the 
surroundings and active malaria drug as this will 
reduce the contact rate between human and mosquito 
through reduction in mosquito population and malaria 
transmission. 

Keyword: Mathematical Model, Next Generation 
Matrix, Sensitivity Analysis, Malaria, Disease Free 
Equilibrium, Basic Reproduction Number. 

I. INTRODUCTION 

Malaria disease is one of the serious and sometimes 
fatal, caused when an infected mosquito carrying the 
parasite feeds on human. Once that infected mosquito 
bites a human, the parasite is transmitted and it 
multiplies in the host liver before infecting and 
destroying their red blood cells. People who are 
infected of malaria are usually very sick and show 
symptoms such as high fever, shaking chills and flu- 
like illness. Anopheles mosquito is the type of 
mosquito that transmits this parasite. Four kinds of 
malaria parasites infect humans: Plasmodium 
falciparum, Plasmodium vivax, Plasmodium ovale 
and Plasmodium malaria. Some of these parasites 


cause severe type of malaria than others. They are 
found in different parts of the world. 

Plasmodium falciparum is the type of malaria that its 
attack usually results in severe infection and if not 
promptly treated may lead to death. In some places, 
malaria can be treated and controlled with early 
diagnosis. However, some countries lack resources to 
do this effectively. 

Prevention of malaria in individuals will generally 
involve the reduction of human mosquito contact 
through the use of bed nets, repellents, prospective 
clothing and mosquito avoidance. This can also 
include reducing the number of mosquito in houses 
through the use of screen and insecticides. Wide scale 
use of mosquito insecticide treated nets has 
contributed substantially to the fall in malaria 
morbidity and mortality. In addition to protecting the 
user, insecticide treated nets protect the community by 
killing anopheles mosquitoes (Moonen et al, 2010; 
Godfray, 2013) and should be deployed in all area 
where malaria is endemic. In as much as this control 
method is very effective, we should observe that the 
main mosquito vectors bite outside in the evening or 
morning and so the protection effect is small. 

Indoor residual spraying with insecticide that persist 
and kill mosquitoes is an important component of 
malaria control (Bernes et al, 2005; Enayati and 
Hemingway, 2010; Pluess et al, 2010). World Health 
Organization (WHO) has developed a global strategy 
frame work for integrated vector management, which 
advocates for an evidence - based, integrated 
approach for vector control and offers guidance for 
successful operationalisation of different approaches 
(WHO, 2004; Beier, 2008). 
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Control activities at community level can utilize 
approaches which directly reduces human mosquito 
contact as well as approaches which reduces the total 
number of mosquitoes in an area. Such approaches 
include the reduction in mosquito breeding ground 
(environmental modification), target the larva stages 
with chemical or biological agents, and massive 
insecticide spraying for the adult mosquitoes. 
Biological control method includes the introduction of 
fish which eat the mosquito larvae or Bercillus 
Thuringiensis which excrete larval toxins. Case 
detection and treatment is another potential control 
method. Identifying and treating infected persons, 
especially asymptomatic individuals, will reduce the 
size of the parasite reservoir within the human 
population and can lower transmission rate. 

II. Related Literatures 

To improve malaria control, there is need to study and 
identify the mechanisms by which other sensory cues 
are detected by mosquito and the potential strategies 
to block them (Agusto et al, 2013; Buonomo, 2013; 
Chitnis, 2010; Hawlay, 2003; WHO, 2013, 2014). 
Insect repellents and other protective odorants can be 
used to prevent bites. Several studies suggested that 
people can acquire protection against mosquito bites 
through regular use of insecticide treated bed net. On 
the other hand, it is challenging to control malaria 
using bed nets alone because bed nets are only used 
for fraction of each day. It is important to examine 
how to combine bed net usage with other control 
approaches affects disease spread. To influence 
mosquito bias requires a combination of disease 
control strategies, such as fumigation of the 
surroundings, repellants, and bed nets, but the 
effectiveness of such a multifaceted approach has not 
been studied using mathematical models. 
Mathematical models for transmission dynamics of 
malaria are useful in providing a better knowledge of 
the disease, to plan for the future and consider 
appropriate control measures. Models have played 
great roles in development of the epidemiology of the 
disease. The study on malaria using mathematical 
model originated from works of Ross 
(1910).According to Ross, if the mosquito population 
can be reduced to below a certain threshold then 
malaria can be eradicated. MacDonald (1957) 
modified the model by including super infection and 
showed that reducing the number of mosquitoes has 
little effect on epidemiology of malaria in areas of 
intense transmission. Dietz and Molinoaux (1974) 
added two classes of humans in their model and these 


includes those with low recovery rate (more 
infectious, greater susceptibility) and those with high 
recovery rate (less infectious, less susceptibility). 
Compartmental models of malaria and differential 
equations are constructed to model the disease 
(Tumwine et al, 2007; Chiyaka et al, 2007; Chitnis et 
al, 2006; Yang, 2000; Ngwa and Shu, 2000). 

Chitnis et al (2006) performed a bifurcation analysis 
of malaria model. Malaria transmission model which 
incorporate immunity in the human population had 
been studied (Tumwine et al, 2007; Chiyaka et al, 
2007; Ngwa and Shu, 2000). Panday et al (2013) 
developed and analyzed a nonlinear SIS mathematical 
model for transmission of malaria caused by infected 
mosquito population on susceptible human population 
and using chemical to control the disease. They found 
that their model exhibits three non-negative equilibria 
which are mosquito and disease free equilibrium E 0 , 
mosquito persistence and disease free equilibrium E 
and endemic equilibrium E*. El-Nor Osman and Adu 
(2017) considered SEIR-SEI model of malaria 
transmission between human and mosquitoes and find 
that malaria may be controlled by reducing the 
contact rate between human and mosquito , the use of 
active malaria drugs, insecticides and mosquito 
treated nets can also help to reduce mosquitoes 
population and malaria transmission respectively. 

Olayini and Obabiyi (2013) used a system of seven 
dimensional ordinary differential equations to study 
the transmission of malaria between humans and 
mosquitoes with nonlinear forces of infection in form 
of saturated incidence rates produce antibodies in 
response to presence of parasite causing malaria in 
both human and mosquito populations. They 
investigated the stability analysis of disease free 
equilibrium and observe that it was asymptotically 
stable when R 0 is less than one and unstable when R 0 
is greater than one. Otieno et al (2016) formulated a 
mathematical model for the dynamics of malaria with 
four - time dependent control measures in Kenya. 
These are insecticide treated bed nets, treatment, 
indoor residual spray and intermittent preventive 
treatment of malaria in pregnancy. They considered 
constant control parameters and calculated the 
reproduction number of disease free and endemic 
equilibrium, followed by the sensitivity analysis and 
numerical simulations of the optimal control problems 
using reasonable parameters. In their study, 
when R 0 < 1, it implies that the disease free 
equilibrium is globally asymptotically stable, but 
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when R 0 > 1, it shows that the endemic equilibrium is 
globally asymptotically stable. Also, at R 0 — 1, their 
model exhibits backward bifurcation and they 
concluded that control programs that follow these 
strategies can effectively reduce the spread of malaria 
disease in different malaria transmission setting in 
Kenya. 

In this paper, we presented the malaria transmission 
dynamics model using ordinary differential equations 
made up of fifteen compartments. Here we considered 
four possibilities in the method which includes the use 
of insecticide treated mosquito nets, treatment, gain of 
immunity and fumigation (indoor residual spraying 
with insecticides, reduction in mosquito breeding 
ground, and targeting the larva stage s with chemical 
or biological agents), after which the model is 
simulate and analyzed. 

III. Model Formulation 

3.1. THE ASSUMPTIONS OF THE MODEL: 

> We assume that treated individuals recover from 
the infection (Chiyaka, 2013); 

> Individuals can acquire protection against 

mosquito bite through regular use of insecticide 
treated net (Agusto et al, 2013; Buonomo, 

2013; Hawlay et al, 2013); 

> Both susceptible, Exposed (Infected), Infectious 
and Recovered humans die at the same natural rate 

> Mosquitoes are infected after biting an infectious 
human (Olaniyi and Obabiyi, 2013); 

> Both susceptible, Exposed (Infected) and 
Infectious mosquitoes die at the same natural rate 
Hh, 

> Infectious human can recover by the effect of the 
natural immune response (Liehl et al, 2015); 

> Treatment may fail due to drug resistance (WHO, 
2010; Phyo et al, 2012). 

3.2. The Parameters are defined as follows: 

p h Recruitment level of susceptible human; 
p m Recruitment level of susceptible mosquito; 
a h Rate at which susceptible human use mosquito 
net; 

a(l — a h ) Rate at which susceptible human 

doesn’t use mosquito net; 

P h Rate at which susceptible human use mosquito 
net and fumigate their surroundings; 
a 1 (1 — [] h ) Rate at which susceptible human who 
use mosquito net doesn’t fumigate their 
Surroundings; 


y h Rate at which susceptible human who doesn’t 
use mosquito net fumigate their surroundings; 
a 2 (l — Yh) Rate of susceptible human who doesn’t 
use mosquito net and doesn’t fumigate 
their Surroundings; 
p h Natural death rate of human; 
p m Natural death rate of mosquito; 
b, b 1 , b 2 Biting rate of mosquito; 
ip lhl ip 2 h> } p 3 h Probability that a bit by an infectious 
mosquito results in transmission of the 
disease to Human; 

xp m Probability that a bit results in transmission of 
the parasite from infectious human to a 
susceptible mosquito; 

Rate at which exposed (infected) human move 
to infectious class; 

Rate at which exposed (infected) mosquito 
move to infectious class; 

k Rate at which exposed (infected) human 

obtain treatment; 

k 1 Treatment failure in exposed (infected) 

human; 

k 2 Treatment failure in infectious human; 
v 1 Rate at which treated exposed (infected) 

human move to susceptible class; 
v 2 Rate of exposed (infected) human who recover 
due to prompt immune response; 
v Rate at which recovered human move to 

susceptible class due to loss of immunity; 
Disease induced death rate of human; 

Disease induced death rate of mosquito; 

8 h Rate at which infectious human are treated; 

(jo h Rate at which treated infectious human move 
to recovered class; 

d Rate at which susceptible human who use 

mosquito net and fumigate their surroundings 
move to susceptible class; 

r Rate recovery of infectious human due to 

immune response; 

p m Rate at which mosquito die due to fumigation; 

The Variables (Compartments) are: 

S h Susceptible human; 

S m Susceptible mosquito; 

S hMN Susceptible human who use mosquito net; 
S h wMN Susceptible human who doesn’t use mosquito 
net; 

S h MNC Susceptible human who use mosquito net and 
fumigate their surroundings; 

ShwMNC Susceptible human who doesn’t use 
mosquito net but fumigate their 
surroundings; 
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S h MNWc Susceptible human who 

use mosquito net 

Em 

Exposed (infected) mosquito; 

but doesn’t fumigate their surroundings; 

h 

Infectious human; 

Skwmnwc Susceptible human 

who doesn’t use 

Im 

Infectious mosquito; 

mosquito net and 

doesn’t fumigate 

I h T 

Infectious but treated human; 

their surroundings; 


R h 

Recovered human; 


E h Exposed (infected) human; 

E h T Exposed (infected) but treated human; 
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We obtain a fifteen dimensional system of ordinary differential equations which describes the progress of the 
disease and its control as: 


dS h 

cit 


Ph F dSfrMNc F ViEhT + vRfr CCh^h Ph$h ^(1 ^h)^h 


dSi 


hMN 


dt 


ttSh. Ph^hMN ^l(l Ph ) $h 


MN 


Ph^hMN 


dS, 


hWMN 

dt 


— a (l a h)Sh. Yh^hWMN a 2 (1 Yh)$hWMN Ph^hWMN 
dShMNC 


dt 


— Ph$hMN ~ dS hMNC Ph^hMNC 


dS , 


hMNWC 

dt 

dShWMNC 


dt 


— a l(l — Ph)ShMN + d 2 I m Pi h^ hMNWC Ph$hMNWC 
~ Yh.Sh.WMN + dim ~ p2hShWMNC ~~ Ph$hWMNC 


dS, 


h WMNWC 

dt 


— a z0- Yh)$hWMN + bl Im Psh^hWMNWC — Ph$hWMNWC 


dE h 

dt 


— Plh^hMNWC + p2hShWMNC + p3hShWMNWC + ^1 EhT ~ kE h — V 2 E h — £ h E h — tl h E h 


dE h T 

dt 


= kE h - v x E h T - \i h E h T - k^bT 


dlu 

+ k 2 I h T ~ rlh ~ %hlh — bin — d h l h — /.i h I h 


dI h T 


dt 
dRh 
dt 
dSrr, 
dt 
dE„ 


dt 


— 8 h lh ~ M h I h T — n h I h T — k 2 l h T 


= ^hPhk + rl h F v 2 E h - vR h - {i h R h 
— Pm F blfi — Ipm^m Pm^m ~~ Pm^m 


Pm^m Cm^rn Pm^m Pm^ri 


dirt 


dt 


(rn^rn b 2 I m b j I m bl m Pmdm Pm Ir 


(3.1) 


From equation (3.1), we observe that the study of the transmission, spread and control of malaria is in two 
interacting population that human which is the host and mosquito which is the vector. The human and mosquito 
population size at the time t are respectively denoted as 

Nh(P) Sfr (t) + SbMN (0 F ShWMN (0 + SfiMNC (0 "F ShWMNePP F ShMNWC (P) F 

S hWMNWC (f) + E h (t) + E h T(t ) + 4(t) + l h T(t) + R h (t ) 

N m (t) = S m (P) + E m (t) + I m (t) 


3.3. Existence and Positivity of Solutions 

In equation (3.1), all the parameters are non-negative and we assume stable population with per capita 
recruitment of susceptible human, susceptible mosquito, death rate of both natural and disease induced and 
those induced on mosquito due to fumigation. For simplicity, the population size of both human and mosquito 
are normalized to 1 each that is S h , Sb MN , S^wmni ShMNC > ShWMNC > ShMNwe » ShWMNWcEh> F/jT, In > IhJ > Rh 
are respectively the fraction of the Susceptible human, Susceptible human who use mosquito net, Susceptible 
human who doesn’t use mosquito net, Susceptible human who use mosquito net and fumigate their 
surroundings, Susceptible human who doesn’t use mosquito net but fumigate their surroundings, Susceptible 
human who use mosquito net but doesn’t fumigate their surroundings, Susceptible human who doesn’t use 
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mosquito net and doesn’t fumigate their surroundings, Exposed (infected) human, Exposed (infected) but 
treated human, Infectious human, Infectious but treated human and Recovered human. Whilc5’ m , E m , I m are 
respectively fraction of Susceptible mosquito, Exposed (infected) mosquito and Infectious mosquito. Here we 
shall show that the system is epidemiologically and mathematically well - posed in the feasible region T given 

b y 


r = H h x H m c R^ 1 x R+ 


Where 


^ and 


Hh — ((S/p ShMN> ShWMN > ShMNC > ShWMNC > $hMNWC > ShWMNWC’ Eh> EhJ > > 4T) ^ ®L+ 1: N h < 

H m = f(S mi E m> l m ) ERl-.N m < (H 

Theorem 1: There exists a domain T in which the solution set 

{$h> ShMN> $hWMN> $hMNC > $hWMNC > $hMNWC > $hWMNWc> EhJ > 61 > IhT> S m , Em> ImJ is contained and 

bounded. 

Proof: Let the initial conditions be 

Shi 0) — Sho,S hMN {0) — S hMNO ,S hWMN (0) — S hWMN()l S hMNC (0) = S hMNCQ , S hWMNC (0) — S hWMNC0 , 

ShMNWci®) — ShMNWCO’ ShWMNWci 0) = ShWMNWCO’ Eh (0) = E h0l E h T(0) — E h T o ,I h (0) — I h0 , 

IhT(Q) — IhJa>Rh(SE) — RhO’ SfniSO = S m 0’ E m i0) = E mO ,I m i0) = Vo (3-2) 

Given the solution set { S h , S hMN , S hWMN , S hMNC , S hWMNC , S hMNWC , S hWMNWC , E h , E h T , I h , I h T , S m , E m , I m \ 
with positive initial conditions (3.2), the human population can be represented as 

Shit) + ShMN (0 + ShWMN (0 + ShMNC (0 + ShWMNC it) + ShMNWC it) + ShWMNWC (0 
+E h (t) + E h T(t) + I h (t) + IhT it) + R h (t) = 1 (3.3) 

Also the mosquito population can be represented thus 

Smit) + Emit) + lmit) = 1 (3.4) 

We can omit the equation for R h in our analysis so that equation (3.1) becomes 
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dS h 

dt 


Ph 3~ ^4mnc 3~ ^1 EfrT <^h^h Ph^h ^/i)4 


+v(l —S h — S hMN — S hWMN S hMNC — S hWMNC — S hMNWC — S hWMNWC — E h — E h T I h — I h T ) 
dS, 


>hMN 

dt 


dSh PhShMN di (1 /4)4mN Ph^hMN 


dS, 


hWMN 

dt 


— d(l — <X h )S h — Yh$hWMN ~ a 2(l — Yh)ShWMN Ph^hWMN 


dS , 


hMNC 


dt 


Ph^hMN dS t 


hMNC 


Ph^hMNC 


dS, 


hMNWC 

dt 

dShWMNC 


— a l(l /4)4 m« + 4 4i x PlhShMNWC Ph^hMNWC 


dt 


— YhShWMN + bl m ^p2h^hWMNC Ph^hWMNC 


dS, 


h WMNWC 

dt 


— a 2 (l — Yh)$hWMN + 44i ~ ^Psh^hWMNWC ~ Ph$hWMNWC 


dE h 

dt 


— 1 Plh^hMNWC + ^2 h4 


hWMNC 


+ ^p 3 h^hWMNWC + 4 EhT ~ kE h — v 2 E h — ( h E h — /. i h E h 


(3.5) 


dE h T 

dt 


= kE h - v t E h T - n h E h T - k t E h T 


dlh. 

dt 


4 Eft T k'2 l h T flh 44 hih 5 h I h Phlh 


dI h T 


dt 
dRh 
dt 
dSrr, 
dt 
dEr, 


dt 


— dhlh ~ d>h!hT — PhlhT — k 2 I h T 


= UhlhT + rlh + V 2 E h - vR h - fi h R h 
— Pm 3" blf j — 1p m ^m ~ Pm^m ~ Pm^m 


(m 4i Pm^m Pm^n 


d I , 


dt 


— 4i4i ~ %mlm ~ b 2 I m — b-^1^ 


bl r 


Pmlm Pmdr, 


Then the time derivatives of H h (t)andH m (t ) along solutions of system (3.1) for humans and mosquitoes 
respectively are denoted by 

4l(0 — Sh(t) + ShMN (0 + ShWMN Id) + SfiMNC CO + OlWMNcCO + S hMNWC (t ) + S hwMNWC (t) 

+E h (t ) + E h T(t ) + 4(0 + 440 (3.6) 

4i(0 — Ph 3” ^4mNC 3” ^1 EhT dfrSfr PhSh — Cl(_l 

+v(l — S h — S hMN — S hWMN — S hMNC — S hWMNC — S hMNWC — S hWMNWC — E h — E h T — I h — l h T) 

+aS h — PhShMN — a l(l — Ph)ShMN — PhShMN 3" a (l — a h)Si l ~ YhShWMN a 2(1 Yh)SfiWMN 
PhShWMN 3" PhShMN ~ dS hMNC — PhShMNC + a l(l — 4)4 mN + b 2 I m ~ l PlhShMNWC ~ PhShMNWC 
3" YhShWMN + bl m ~ U*2hShWMNC PhShWMNC + a 2 (1 — Yh)ShWMN + b 1 I m ~ x l J 3hShWMNWC 
PhShWMNWC 3" 'PlhShMNWC 3" 4 , 2hShWMNC 3" 'P'ihShWMNWC 3" k-iE h T — kE h — V 2 E h — <441 — Ph^h 

+kE h — fi h I h T — v t E h T — k t E h T + <44. + b 2 I h T — n h I h — rl h — 44 — ^4 — 44 + Shlh — <y /i47’ 

~PhIhT ~ k 2 IhT 


+ Sh 


It is worthy to note that in the absence of the diseas ebl m , b-^Ifn, b 2 I m , v 2 E h ,rI h , 44- bI h anda) h I h T will all be 
zero. Therefore we have 
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Hfi (0 — Ph — iph^h 4 " Ph^hMN 4 " Ph^hWMN 4 " Ph^hMNC 4 ~ Ph^hWMNC 4 " Ph^hMNWC 4 " Ph^hWMNWC 4 ~ Ph^h 

+ PhEh T + A0i4 + PhlhT) + ^(1 - Hh) 

H h (t) = p h - p h H h + v(l - H h ) 

H h (0 = Ph + v - Qi h + v)H h 

H h (t) + (p h + v)H h <p h + v (3.7) 

Integrating both sides of equation (3.7) using integrating factor method (Kar and Jana 2013; Birkhoff and 
Roffa, 1989; Chinebu et al, 2018) we have 

H h +p(t)dt = f{t) 


Where 


H h < e -/p(0dt (J e fpW dt f(t)dt + c) 
Pit) = p h + vandfit ) = p h + v 


We then have the integrating factor as 

r(t) = e/pOJdt = e f(nh+v)dt _ e (nh+v)t 

Now integrating (3.7) using the value of the integrating factor r(t) = we have 

1 


H h it) < 


r(t) 


I 


r(t). f(t)dt + c 


H h it ) < 

< 


e (n h +v)t 

1 


I 


,(n h +v)t 


. (ph + v)dt -I- c 


e (Ph+v)t 

1 

< 


ip h + v) j e^ h+v ^ t . dt + c 
iph + v) 


(, Hh+v)t 


H h (t ) <^i— + ce -0th+»)t 
Ph + v 


iPh + v) 


,(Hh+v)t + 


(3.8) 


Where c is the constant of integration and if we let t -» oo we have 

Ph + v 


But 


Eh ” 1 ” $hMN 4 ” $hWMN 4 ” $hMNC 

PhiPh + d 2 + Ph)id + Ph) 


H h it) = 

Ph + v 

4" ShWMNC 4” ShMNwe 4" ^hWMNwe 4” Eh 4- E h T + 7/j + /^T 


S h < 


h iPh + a h 4- djifih 4" d 2 4- p h )id 4- Ph) ~ dp h tth 
which at disease free equilibrium gives 


(3.9) 


S h = 


Ph 


Also, 


a + Ph 

Hmit) = Emit ) + Emit) + 7 m (t) 

H-mit) — Pm 4” dlh — Pm$m ~ Pm$m 4” fm^m (Vn^Vn — Pm^m — PmE-n 


+(mE m ~ %mlm ” b 2 I m ~ b i/„ 


bL 


Pmdm Pm P, 


(3.10) 


But in the absence of the disease, b 2 I m , h l I m , bl m , P m I m , p m S m , p m E m , p m I m will all be zero and we obtain 

Hmit) — Pm — Pmi^m 4” E m + I m ) 

If we let S m + E m + I m — H m , then we have 


Where 


Hmit) Pm PhHm 

Hmit) + PmH m < Ph (3.11) 

H m + pit)dt = fit ) 

H m < e~SPd)dt ^ j efpWdtf^fo + 

Pit) = p m andfit) = p m 
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We then have the integrating factor as 

_ £>/ p(t)dt g/ Ll m dt gtPmP 

Now integrating (3.11) using the value of the integrating factor r(t) = we have 

H m {t) < -!- [ f r(t).f(t)dt + c x 
r(t) U 


Hmit) < 


> (Mm)^ 


/ 


i (M??i)^ 


■(Pm)dt + Ci 


< 


£ (Mm)£ 
< 


> (Mm)^ 


(Pm) | e^ m)t .dt + c ± 
PjH}L P (Pm)t 


H m (t ) < — + c x e 
Pm 

Where c x is the constant of integration and if we let t -> oo we have 


Pm 


+ Ci 


(3.12) 


But 


Pm 

limit) ~ ~~ = S m + E m + I m 

r-m 


Pm 


which reduces to — at disease free equilibrium. We observe from the dynamics described by thee system (3.1), 

Pm 

(3.9) and (3.13) that the region 

r = | iSh> $hMN> SfiWMN> ShMNC > $hWMNC > $hMNWC > hWMNWC> ^h> E'hJ > Ifi > IhT< $m> E m > lm) £ ^+ 4: 


Pm ~b Tm "b Pm 


(3.13) 




Um ) 


& "b Ph Pm J 

is positive invariant. This simply implied that the system for humans and mosquitoes are respectively well 
posed mathematically. Then for the initial stating point H h c M+ 1 and H m c M+, the trajectory lies in T. Thus, 
we focus our attention only on the region T. 


3.4. Stability of the Disease Free Equilibrium (DFE) 

The disease free equilibrium state denotes the point at which the differential equations of system (3.1) are equal 
to zero and it refers to absence of infection. Therefore, all the infected class will be zero and the entire 
population for human will comprise of Susceptible human, Susceptible human who use mosquito net, 
Susceptible human who doesn’t use mosquito net, Susceptible human who use mosquito net and fumigate their 
surroundings, Susceptible human who doesn’t use mosquito net but fumigate their surroundings, Susceptible 
human who use mosquito net but doesn’t fumigate their surroundings, Susceptible human who doesn’t use 
mosquito net and doesn’t fumigate their surroundings while the entire population for mosquito will be the 
susceptible mosquito only. 


At equilibrium state we set the rate of change of each variable to be equal to zero, that is 
dS h dS hMN dS hWMN dS hMNC dS hMNWC dS hWMNC dS hWMNWC dE h 


dt 


dt 


dlfi 

dt 


dt 

_ dI h T 

dt 


dt 

dR h 

dt 


dt 


dt 


dt 


= 0 


dt dt 
(3.14) 


dE h T _ 

dt 


Let (S h , S hMN , S hWMN , S hMNC , S hWMNC , S hMNWC , S hWMNWC) E h , E h T , I h , I h T, S m , E m , I m ) 

_ / f 1 ^ P ^ ^ ri ^ ri t t n ♦ t~» H* rr 1^1 P ^ T~> ^ t 

— Wh >^hMN > b hWMN >t>hMNC > b hWMNC hMNWC > b h.WMNWC > > ^h 1 > l h > l h l > b m < h m > l m ) 

at equilibrium state. If we substitute E h * — E h T* — I h * — I h T* — E m * — I m * — 0 into equation (3.1) we have 
the disease free equilibrium state for human and mosquito calculated as 
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Ph 4” dShMNC 4” ViEfrT + vRfr — oc^S^ — Ph$h — a(l — 0Ch)Sh — 0 ^ 

a $h — Ph^hMN — a l (1 — Ph)ShMN ~ Ph$hMN ~ 0 
a (l — a h)Sh ~ Yh.Sh.WMN <^2(1 Yh)$hWMN ~ Ph^hWMN ~ 0 
Ph$hMN ~ dS hMNC — PhShMNC ~ 0 
a l(l — Ph)ShMN 4" d 2 I m ~ 'PlhShMNWC ~ Ph^hMNWC ~ 0 
Yh$hWMN 4" him ~ ^p2hShWMNC ~ Ph^hWMNC ~ 0 
a 2(l — Yh)$hWMN 4" bllm ~ x l } 3h$hWMNWC ~ Ph$hWMNWC ~ 0 

1 Plh$hMNWC 4" ^p2hShWMNC + x P'3hShWMNWC 4" k x E h T — kE h — V 2 E h — Y, h E h — p h E h — 0 (3.15) 

kE h - v x E h T - p h E h T - k x E h T = 0 
(h^h + k 2 I h T — rl h — % h I h — bl h — S h I h — p h I h — 0 
f/i — <^/i//iT — Ph^hT ~ k 2 IhT — 0 
a) h I h T 4- rl h + v 2 E h - vR h - p h R h = 0 
Pm 4” 5//! — — Pm^m ~ Pm^m ~ 0 

^Pm^m Cm Em ~ p m E m — /i m E m — 0 

Cm^rn — Cm^m — ^2 I m — b x I m — bl m — p m I m p m I m — 0 j 

If we let a( 1 — a h ) — d 1 a 1 ( 1 — (3 h ) — d 2 a 2 (l — y h ) = d 3 and label equation (3.15) as itoxv, we calculate 
the equilibrium point for both human and mosquito. 

Substituting E h * — E h T* — I h * — I h T* — E m * — I m * — 0 we have from (3.15xii) as 

u h IhT 4- rl h + v 2 E h - vR h - p h R h = 0 

But I h T — I h — E h — 0 

(y + p h )R h = 0 =>R h = 0 

From (3.15iv) we have 


(3.16) 


>hMNC 


Ph$hMN ~ dS hMNC 
Ph$hMN 

d + p h 


~ Ph^hMNC — 0 


(3.17) 


Equation (3.15ii) gives 


Ph + d 2 4 - p h 
Substituting (3.18) into (3.17) we obtain 

c _ Ph^h^h 

^hMNC 


a h$h ~ Ph$hMN ~ a l (1 — Ph)$hMN ~~ Ph$hMN ~ 0 
OChSh 

ShMN = W , i , .. (3.18) 


(d + Ph)(Ph + d 2 + p h ) 
Substituting (3.16) and (3.19) into (3.15i) gives 

dPh^h^h 

Ph + 


(3.19) 


( a h 4 - d x + Pfi)S h - 0 


(d + Ph)(Ph + d 2 + p h ) 

Ph(d + PfJiPh + d 2 4- p h ) -1- [ d/3 h a h - (a + d x + p h )(d + p h )(P h 4- d 2 + p h )]S h = 0 


Substituting (3.20) into (3.18) gives 


a h Ph(.d+p h ) 


hMN ( a h + d x 4 - p h )(d + p h )(P h + d 2 + p h ) - dp h a h 

From (3.15iii) we have 

a (l — a h)Sh ~ Yh$hWMN a 2 (1 Yh)ShWMN ~ Ph^hWMN ~ 0 

diP h (d + p h XP h + d 2 + p h ) 


(3.21) 


>hWMN 


Substituting (3.21) into (3.17) gives 

ShMNC ~ 


(Yh + d 3 + p h )[(a ,i + d x + p h )(d + p h )(Ph + d 2 + p h ) - dp h a h ] 

Ph^hPh 


(3.22) 


(a h + d x + Ph)(d + p h XPh + d 2 + p h ) - dp h a h 


(3.23) 
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From equation (3.15v) we have 

a l(l — Ph)ShMN + ^2An — H’lh^hMNWC ~ Ph^hMNWC ~ 0 
d 2 cf/iP/i(d + pX) 

Jh 


>hMNWC ~ 


(3.24) 


0 pih + Ph)K«h + rfi + fihXd + PhXPh + d 2 + Ph) - dp h a h ] 

From equation (3.15vi) we have 

Yh.Sh.WMN + bl m ~ IpZhShWMNC ~ Ph$hWMNC ~ 0 

= _ YhdiPhid + p, h XPn + d 2 + p h ) _ 

hwMNC ^ + M/i )[( 7/i + d 3 + j u /l )[(a /l + d 1 + p h )(d + p h )(J3 h + d 2 + Hh) ~ dp h a h ]] 
From equation (3.15vii) we have 

a 2 (l — Yh)ShWMN + bil m — IpshShWMNWC ~ Ph$hWMNWC ~ 0 

= _ d 3 d 1 p h (d + Ph)(Ph + d 2 + p h ) _ 

hWMNwc ^ M h )[(y h + d 3 + Ph)K a h + rfi + Ph)(d + p h )(J3 h + d 2 + p h ) - d/S h a h ]] 
From equation (15xiii) we have 


(3.25) 


(3.26) 


c _ P m 

‘-’m 

Pm 


Pm 3” b^h ^Pm^m Pm$m Pm$m 0 


(3.27) 


3.5. Analysis of the System of Control 

We analyzed the equilibrium point and the basic reproduction number of the systems having in mind that the 
control parameters a h (t),(3 h (t),y h (t),8 h (t),k(t) and p(t) are constant functions. 

3.6. Equilibrium Point and Basic Reproduction Number R 0 

Here we calculate the disease free equilibrium point of system (3.1) which shows that the population is free of 
malaria infection and this implies that there isn’t need for control. Therefore, the control parameters will be 
equal to zero. 

3.7. Disease Free Equilibrium 

Let M 0 * be the disease free equilibrium point for both human and mosquito, then we have 


M 0 — (ShO’SfiMNO’Sh 


hWMNO’ ShMNCO’ J hWMNCO> J hMNWCO< °hWMNWC0 


)> s. 


,s h 




3,0,0, 0, 0,0, S m o, 0,0 ) 


Remember that 


and 


S h0 — 


Ph(d + Ph)(Ph + d 2 + p h ) 


(a + d 1 + p h )(d + p h )(J3 h + d 2 + p h ) - d/3 h a h 


a( 1 - a h ) = d x ,%(1 - p h ) = d 2 a 2 (l - y h ) = d 3; 

c _ Ph 

-j ho — ; 

a -I - Ph 

.. __ ahPh(d + Ph) _ 

hMNo ^ a ^ + d i+ p h )(d + p h )((3 h + d 2 + p h ) - dp> h a h 

dxPnid + p h XPh + d 2 + p h ) 


= 0 


5 hWMNO ~ 


C Yh + d 3 + P h )[(a h + dj + Ph)(d + p h )Q3 h + d 2 + p h ) - d/3 h a h ] 
c a Ph 


’hWMNO 


(a + p h )(a 2 + p h ) 

Ph^-hPh 


>hMNC0 


(a h + d 1 + p h )(d + p h X(3 h + d 2 + p h ) - d(l h a h 
YhdiPhid + p h XP h + d 2 + p h ) 


= 0 


hwMNC 0 (xp 2h + p h )[(y h + d 3 + p h )\(a h + d x + p h )(d + p h )(J3 h + d 2 + p h ) - d(l h a h ]] 

d 2 a h p h (d + Ph ) 


= 0 


5 hMNWCO 


= 0 


’hWMNWC 0 


0 Pih + PhX&h + di + Pn)(d + p h )(/3 h + d 2 + p h ) - d/3 h a h ] 
_d 3 djp h (d + p h )(p h + d 2 + p h ) _ 

0/>3/, + Ph)[(.Yh + d 3 + p h )[(a h + d x + p h )(d + p h )(J3 h + d 2 + p h ) - dp h a h ]] 
_ aa 2Ph 

Ph( a + Ph)(a 2 + p h ) 
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Pm Pm 
Pm 


C — C 
J m 0 — ~ 


Pm 4" t pm 

The disease free equilibrium point for human and mosquito is 

M a ’ = ( ———,0,7- ^--. 0 , 0 , 0 . aa2Pk 

\a ' ' 


+ Ph’ ’ ( a + Ph)( a 2 + Ph) 


' Ph(a + p h )(a 2 + n h ) 


,0,0,0,0,0, — ,0,o) (3.28) 

Pm ' 


3.8. Basic Reproduction Number 

The basic Reproduction number of epidemiological model usually denoted by R 0 is the average number of 
secondary infections caused by a single infectious individual during their entire infectious lifetime (Driessche 
and Watmough, 2002). It helps to show whether the infection will be endemic or not i.e., spread through the 
population or not. When the basic reproduction number is less than one ( R 0 < 1), the disease free equilibrium 
point is locally asymptotically stable while when it is greater than one ( R 0 > 1), the disease free equilibrium 
point is unstable.To obtain R 0 for model system (1) we use the next generation matrix technique (Van den 
Driessche and Watmough, 2002; Diekmann et al 1990). We shall start with those equations of the model that 
describe the production of new infections and change in state among infected humans and mosquitoes. 

Let Q — [S h , S hMN , S hWMN , S hMNC , S hWMNC , S hMNWC , S hWMNWC , E h , E h T , l h , I h T , S m , E m , I m ] T where T 
denote transpose. Then 

dQ 

FdQ)-V(Q) 


dt 


where 


no) 


0 


a(l - a h )S h 
0 
0 
0 

a z0- ~ Yh)ShWMN 

U’lh^hMNWC + 4*2 h^hWMNC + 4 J 3h^hWMNWC + ^E^T 

kE h 

ShE h + k 2 I h T 

Shh 

4m Em 


and 


-V(Q) = 


(Vl "b Ph)^hMN 
( a 2 + Ph) ShWMN 

(d + Ph)S hM NC PhShMN 
(<Pl h + Ph)$hMNWC ~ a l(l — Ph)ShMN 
(PPzh + Ph)$hWMNC ~ YShWMN 
(<P3 h + Ph) ShWMNWC 
(k + V 2 + p h + 4h)Eh 
Or + k x + p h )E h T 
( r + Ph + b + 4h + 8 h )I h 
(<*>h + ^2 + Ph)Ih T 
(Pm ”b Pm T 4m) 

(b 2 + b t + b + + p m + li m ) 
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If we find the derivatives of F and V at the disease free equilibrium point M 0 we obtain F and V respectively 
as; 


F := 


0 0 
0 0 
0 0 
0 0 


0 


0 


0 0 
0 0 
0 0 
0 0 

0 0 


0 a 2 (l-T h ) 0 0 


0 0 0 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 0 0 0 b„ 


0 0 0 0 0 0 b 


0 0 
0 0 

0 0 

0 0 
0 0 


< ’h 


0 0 0 0 0 0 b 


0 V lh V 2h V 3h k l ° 0 0 0 0 


k 0 0 0 0 0 


0 k 2 0 0 0 


0 0 0 5 U 0 0 0 

n 

0 0 0 0 0 0 0 

0 0 0 0 0 c o 

Tn 


V := 


V^h 

0 

-Ph 

-IK 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

-?h 

0 

0 

0 

0 

0 

0 

0 


0 

0 

d + ^ 


0 


0 


0 


V v ih 

0 vv 2h 


V V 3h 

0 

0 


0 

0 

0 

0 

0 

0 

k + Wti 
0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

V k i + % 
0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

V r + 5 h^h 
0 

0 

0 


H h + C o h + k 2 


H + Q + p 
m tti r n 


|4 + 5 + P 
m m r mJ 
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a , + ^ 


a 2 + h 


i + K +d Wh 




all - 


W + w a iWh 




1 p h 


0 0 


0 0 


0 0 


0 0 


0 0 


K + a 2'V a 2Wh 




0 


k + W^ 
0 0 


k ify“h 


0 0 


F-V 


-1 


0 

0 


0 

0 0 

0 

0 

0 

1 

0 

0 

0 


r+ ^h +ll h +8 h 

0 

0 

0 

0 

0 0 

0 

0 

0 

1 

0 

0 


k 2 + V a h 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

1 

0 


p + L+P 

m m m 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 

1 


l +LI +p 

7n r ra r m 

0 


0 


0 0 

0 

0 

0 

0 0 

0 

0 

0 

0 


0 


0 0 

0 

0 

0 

0 0 

0 

0 

0 


0 0 0 0 


0 0 0 0 


0 0 0 0 0 

0 0 0 0 0 


l + U +P 
’m r m r m 


a (‘-Phh 


v'h 

a 2 + 0 h 


V h 


0 0 0 0 0 0 0 0 0 


0 0 0 0 0 0 0 0 0 


t, + p +p 
’m m r m 


l + \i + p 
r m r m 


0 


Vi v 2 v 3 


k 


M + a iV a iWh w +a 2'V a 2Wh 

0 0 0 0 0 0 


v 1 + p h V% V o h k + v 2 + p h+ i; h 
k 


0 0 0 0 


k+ w^ 


0 0 0 0 


k„ 


0 0 0 0 

0 0 0 0 

0 0 0 0 


0 


k+ W^ r + ^h + V 5 h 

8 h 

0 0 


0 0 


0 0 


r+ w\ 

0 0 0 0 0 


0 0 0 0 0 0 0 0 


o +L + P 

m in r m 
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Observe from determinant of matrices that 

> When every element in a row (or column) of an m x m (square) matrix is zero, the value of the determinant 
will be zero. 


Since the matrix above satisfies this property, we state that the 

|F.R _1 | = 0 

To obtain the basic reproduction number R 0 we use the trace and determinant of the matrix and we let FV~ 1 — 
A. 


Let 


1 

R 0 — <p(i4) = - {trace (A) + Jtrace (A) 2 — 4det(A)) 


(3.29) 


Since det(A) = 0, we have that 


R 0 


k + v 2 + p h + Zh r + Pi, + Kh 


(3.30) 


Therefore, if R 0 < 1, then the disease free equilibrium is stable and malaria infection can be controlled using 
the three control measure. On the other hand, if R 0 > 1, then the disease free equilibrium is unstable and the 
malaria infection will persists’ in the population. 

3.9. Sensitivity Analysis of the Basic Reproduction Number R 0 of Malaria with Control 

We conducted sensitivity analysis to determine the nature of the model. R 0 is the expected number of new 
infected humans as a result of introducing one infected human in a population at a disease free equilibrium in 
the absence of disease control methods. If R 0 < 1, each infected human infect on average less than one 
individual in the population and hence the disease dies out. But if R 0 > 1, each infected human infects more 
than one individual in the population and thus the disease is able to invade the susceptible population. In other 
words, R 0 is used here to determine the effectiveness of the control measure. To calculate the value of R 0 we 
used the parameter values in table 1. 


Table 1 Parameter values for the calculation of R 0 and Numerical Simulation of the Model. 


Parameters 


Value 

Range 

Ph 

Recruitment level of susceptible human 

0.0814 

Pm. 

Recruitment level of susceptible mosquito 

0.0184 

a h 

Rate at which susceptible human use mosquito net 

0.75 

a( 1 - a h ) 

Rate at which susceptible human doesn’t use mosquito net 

0.25 

Ph 

Rate at which susceptible human use mosquito net and fumigate their 

surroundings 

0.1125 

Yh 

Rate at which susceptible human who doesn’t use mosquito net fumigate their 

surroundings 

0.11 

Ph. 

Natural death rate of human 

0.0000548 

Pm 

Natural death rate of mosquito 

0.67 

b 



bhl,2,3 

Biting rate of mosquito 

0.12 

^Phl,2,3 

Probability that a bit by an infectious mosquito results in transmission of the 

disease to Human 

0.0200 

l 'Pm 

Probability that a bit results in transmission of the parasite from infectious human 

to a susceptible mosquito 

0.5 or 
0.8333 

(h 

Rate at which exposed (infected) human move to infectious class 

0.0588 

("m 

Rate at which exposed (infected) mosquito move to infectious class 

0.0017 

k 

Rate at which exposed (infected) human obtain treatment 

0.28 
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fcl 

Treatment failure in exposed (infected) human 

0.09 

^2 

Treatment failure in infectious human 

0.017 

Vi 

Rate at which treated exposed (infected) human move to susceptible class 

0.05 

V 2 

Rate of exposed (infected) human who recover due to prompt immune response 

0.00137 

V 

Rate at which recovered human move to susceptible class due to loss of immunity 

0.000137 

U 

Disease induced death rate of human 

0.02326 

%m 

Disease induced death rate of mosquito 

0.07 

8 h 

Rate at which infectious human are treated 

0.3325 

Uh 

Rate at which treated infectious human move to recovered class 

0.000548 

d 

Rate at which susceptible human who use mosquito net and fumigate their 
surroundings move to susceptible class. 

0.1 

r 

Rate recovery of infectious human due to immune response 

0.0035 

Pm. 

Rate at which mosquito die due to fumigation 

0.15 

%(1 -Ph) 

Rate at which susceptible human who use mosquito net doesn’t fumigate their 

surroundings 

0.8875 

0-2(1-Yh) 

Rate of susceptible human who doesn’t use mosquito net and doesn’t fumigate 

their surroundings 

0.89 


To calculate the value of R 0 , we use the parameters as stated in table. 

/Ci ky 

R — _ t _|_I_ 

0 k + v 2 + Hh + Zh r + Ph + Kh + &h 

0.09 _ 0.0175 

R ° ~ 0.28 + 0.00137 + 0.0000548 + 0.0588 + 0.0035 + 0.0000548 + 0.02326 + 0.3325 

0.09 i 0.0175 

” 0.3402248 + 0.3593148 
= 0.2645309807 + 0.0487038107 
R 0 = 0.3132347914 < 1 

> For the effect of /q on R 0 we have 

B dR Q R 0 

E*° = —- x = 0.7280723238 

kl dk x k x 

> For the effect of k on R 0 we have 

dR 0 R 0 

E*° = —- x -f- = -0.8698062743 

u K K 

> For the effect of v 2 on R 0 we have 

d dR 0 R 0 

= —— x — = -177.77061763 

2 OV 2 V 2 

> For the effect of ^ on R 0 we have 

R dR 0 R 0 

E^ 0 = ——- x —- = -4.1419346281 

<A d( h 4 

> For the effect of r on R 0 we have 

□ dR 0 R 0 

E^° = —- x — = -12.130810762 
or r 

> For the effect of k 2 on R 0 we have 

R dR 0 R 0 

E? 0 = —- x = 49.814621989 

k2 dk 2 k 2 

> For the effect of fi h on R 0 we have 

dR 0 R 0 

Eu? = -r— x — = -5219.0436284 
dfi h n h 

> For the effect of <f/, on R 0 we have 
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R dR 0 R 0 

E^ 0 = —- x -E- = -1.8253583204 

^ df h 

> For the effect of 8 h on R 0 we have 

E^° = ^ x ^ = -0.1276927339 
d8 h 8 h 


Table 2: The Effect of the Parameters on R 0 


Parameters 

Effect on R 0 

ki 

0.7280723238 

k 

-0.8698062743 

v 2 

-177.77061763 

4 

-4.1419346281 

r 

-12.130810762 

k-2 

49.814621989 

Eh 

-5219.0436284 

Zh 

-1.8253583204 

8h 

-0.1276927339 


The sensitivity of E(k_l ) and E(k_2 ) are positive and this implies that the value of R_0 increases as the values 
of k_l and k_(2 ) increases. Contrary to this, the values of the parameters 

E(k),E(5_h),E(^_h),E(p_h),E(r),E(^_h) and E(v) are negative, indicating that the value of R_0 decreases as 
k,8_h,2_h,p_h,r,c_h and v increases. This simply shows that these parameters contribute in reducing that rate 
of transmission of malaria infection over time in the population. If these parameters are maintained, the 
transmission of the disease may decrease causing the cases of malaria infection in the population to reduce and 
go below the endemicity threshold. 


IV. Results and Discussion 

In this section, we study numerically the behavior of system (3.1) using the parameter values given in table 1 
and by considering initial conditions, II = (4,(0), S hMN (0), S hWMN (0), S hMNC (0), S hWMNC (0), 

ShMNWc (0), S hWMNWC (0), E h ( 0), E h T( 0), 4(0)47(0), R h ( 0), S m (0), E m ( 0), I m { 0)}. The numerical simulation 
are conducted using Matlab software and the results are given in figure 2-4 where figure 2a - 2b illustrate the 
behavior of the reproductive number for different values of the model parameter k and r respectively. Figures 
3i - 3xv and 4i - 4xv shows the varying effects of when there is no control and when there is control 
respectively. 


The basic reproduction number of the system is given by 

/Ci kn 

+ 


R 0 = 


= 0.3132347914 < 1 


k + v 2 + Eh + 4 r + Eh + 4 + <4 
This simply shows that the basic reproduction number is less than one, hence the disease free equilibrium is 
stable which indicates that the malaria infection can be controlled in the population using the control measures 
stated in the model. However, it also confirms the result of the sensitivity analysis of R 0 in table 2 and we then 
state that with effective use of insecticide treated bed nets, modification of the environment through fumigation 
and treatment of infectious human, the future number of malaria infection cases will reduce in the population. 



k=Control fig2a 


@ IJTSRD I AvailableOnline@www.ijtsrd.coml Volume-2 I Issue - 6 ISep-Oct2018 


Page: 355 






















International Journal of Trend in Scientific Research and Development (IJTSRD) ISSN: 2456-6470 



vh2=Control fig2b 


Figures (2a - 2b): Numerical Simulation of the basic Reproduction Number R 0 using 

different rate of k and r. 
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Figures (3i-3xv): Numerical Simulation of model system (3.1) when there is no 
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Figures (4i-4xv): Numerical Simulation of model system (3.1) when there is 
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From figures (2a - 2b), we observe that the basic 
reproduction number which determines the dynamics 
of the disease in response to the control parameters is 
less than unity. This means that increase in 
k and r can slow disease spread. Eradication is easily 
done if k and r with the other control parameters is 
large relative to k 1 and k 2 . Therefore, increasing 
k and r means increasing the effectiveness of the 
control measures, thus, malaria disease spread can be 
stopped. 

Figures (3i - 3xv) shows the behavior of Susceptible 
human, Susceptible human who use mosquito net, 
Susceptible human who doesn’t use mosquito net, 
Susceptible human who use mosquito net and 
fumigate their surroundings, Susceptible human who 
doesn’t use mosquito net but fumigate their 
surroundings, Susceptible human who use mosquito 
net but doesn’t fumigate their surroundings, 
Susceptible human who doesn’t use mosquito net and 
doesn’t fumigate their surroundings, Exposed 
(infected) human, Exposed (infected) but treated 
human, Infectious human, Infectious but treated 
human and Recovered human together with 
Susceptible mosquito, Exposed (infected) mosquito 
and Infectious mosquito when there is no control. We 
observe that the susceptible human population drops 
because the population of exposed (infected) and 
infectious human increases as the population of 
exposed (infected) and infectious mosquitoes 
increases. 

Similarly, figures (4i - 4xv) depict the varying effects 
of the proportion of the control measures against 
mosquitoes and malaria infection. As shown in figures 
(4i - 4xv), the population of susceptible mosquitoes 
decrease with time as a result of the fact that the 
density of human who use insecticide treated bed nets 
and fumigate their surroundings increases. Also 
observed here is the increase in the density of the 
exposed but treated and infectious but treated humans 
and this indicates that the treatment is effective. We 
therefore state that when sufficient human population 
effectively use insecticide treated bed nets and 
adequately fumigate their environment, his will 
reduce the density of the exposed (infected) humans. 
More so, the application of an effective therapy will 
reduce the density of the infectious humans as 
indicated by the simulation. 


In conclusion, the numerical simulation result showed 
that combination of the use of insecticide treated bed 
nets, fumigation of the surroundings and adequate 
treatment can effectively control malaria infection. 

V. Concluding Remarks 

This paper designed and analyzed a compartmental 
model for malaria transmission dynamics in human 
and mosquito population with combined control 
methods. The human population was divided into 
twelve compartments of Susceptible human, 
Susceptible human who use mosquito net, Susceptible 
human who doesn’t use mosquito net, Susceptible 
human who use mosquito net and fumigate their 
surroundings, Susceptible human who doesn’t use 
mosquito net but fumigate their surroundings, 
Susceptible human who use mosquito net but doesn’t 
fumigate their surroundings, Susceptible human who 
doesn’t use mosquito net and doesn’t fumigate their 
surroundings, Exposed (infected) human, Exposed 
(infected) but treated human, Infectious human, 
Infectious but treated human and Recovered human 
while the mosquito population was divided into three 
compartments of Susceptible mosquito, Exposed 
(infected) mosquito and Infectious mosquito, since 
infectious mosquito does not recover from the 
infection. 

We showed that the system is epidemiologically and 
mathematically well - posed by determining the 
positive invariant and the disease free equilibrium 
point was also shown. The basic reproduction number 
was obtained using the next generation matrix 
technique. We prove that the disease free equilibrium 
is asymptotically stable when the basic reproduction 
number is less than one and unstable when is greater 
than one. 

Finally we numerically simulate our model to see the 
effect of the control parameters. Our result showed 
that increasing the proportion of humans who use 
insecticide treated bed net and fumigate their 
surroundings will go a long way in reducing human - 
mosquito contact rate, thereby reducing the density of 
the exposed (infected) human in the population. 
Similarly, effective treatment of infectious human 
helps to boost their immunity as this will prevent the 
occurrence of malaria infection over time. Therefore, 
combined control is more effective than when there is 
a single control method or when there no control. 
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